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In this paper we extend the range of consistency of a constant bulk viscosity model to redshifts up 
to z ~ 8.1. In this model the dark sector of the cosmic substratum is a viscous fluid with pressure 
p = —^9, where is the fluid-expansion scalar and ^ is the coefficient of bulk viscosity. Using 
the sample of 59 high-redshift GRBs reported by Wei (2010), we calibrate GRBs at low redshifts 
with the Union 2 sample of SNe la, avoiding then the circularity problem. Testing the constant 
bulk viscosity model with GRBs we found the best fit for the viscosity parameter (^ in the range 
< C < 3, being so consistent with previous probes; we also determined the deceleration parameter 
go and the redshift of transition to accelerated expansion. Besides we present an updated analysis 
of the model with CMB5-year data and CMB7-year data, as well as with the baryon acoustic peak 
BAO. From the statistics with CMB it turns out that the model does not describe in a feasible way 
the far far epoch of recombination of the universe, but is in very good concordance for epochs as far 
as 2: ~ 8.1 till present. 
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I. INTRODUCTION 

Since 1998 reliable cosmological data have been ac- 
cumulated leading to the conclusion that our universe 
has entered recently into an accelerated expansion epoch. 
The luminosity-distance data of supernovae la (SNe la) 
has provided the main evidence. Hitherto explanations 
to this phenomenon are not satisfactory. Most models 
consider accelerated expansion as due to a component 
of the universe that behaves opposite to gravity, the so 
called dark energy (DE) . Another unknown component of 
the universe is the dark matter (DM), the missing mass 
necessary to held together galaxy clusters, also needed to 
explain the current large scale structure of the universe. 
The preferred model to describe accelerated expansion is 
the Lambda-Cold-Dark-Matter (ACDM). Hovifever, cos- 
mological constant has several theoretical drawbacks, like 
the discrepancies between the observed and theoretical 
values, i.e. the fine tuning problem; the coincidence prob- 
lem is another objection, meaning the fact that dark 
energy is nearly equal to the matter density just now. 
Among many proposals to give a satisfactory answer 
there are some models that consider DE as a dynamical 
component, such as Quintessence or K-essence; in others 
proposals DE appears as the result of the interaction of 
fundamental particles. In here we adhere to the stream 
that considers DM and DE as different manifestations of 
the same component, describing the dark sector as some 
kind of matter whose physical properties depend on the 
scale: behaving as DM at high densities and transforming 
into DE at lower densities. These Unified Dark Matter 
models have the generalized Chaplygin gas as its paradig- 
matic example [iHSIj however these models present un- 
desirable features like oscillations in the matter power 
spectrum and to mend this inconvenience it should be 
assumed the existence of nonadiabatic pressure pertur- 
bations. 

On the other side, in the context of inflation of the very 
early Universe, it has been known since long time ago 
that an imperfect fluid with bulk viscosity can produce 



an acceleration in the expansion without falling back on 
a cosmological constant or some inflationary scalar field 
|^-|6[. Inspired in those inflationary models there are 
some recent [THlSj and not so recent IJI developments 
that assume a universe filled with a viscous single fluid. 
The bulk viscosity contributes to the cosmic pressure and 
drives the accelerated expansion. 

The origin of the bulk viscosity in a physical system 
is due to its deviations from the local thermodynamic 
equilibrium [1^, [l^ and different cooling rates [I3, Il3l ■ 
In a cosmological fluid, the bulk viscosity arises when 
the fluid expands (or contracts) too fast so that the sys- 
tem does not have enough time to restore the local ther- 
modynamic equilibrium and then it appears an effective 
pressure restoring the system to its local thermodynamic 
equilibrium. Bulk viscosity can be seen as a measurement 
of this effective pressure. When the fluid reaches again 
the thermal equilibrium then the bulk viscosity pressure 
ceases [l3|, [2(]| - |22| . So, in an accelerated expanding uni- 
verse it may be natural to suppose that the expansion 
process is actually a collection of states out of thermal 
equilibrium in a small fraction of time that being mod- 
eled by a bulk viscous fluid might be a more realistic 
description of the accelerated universe today. 

Moreover, these models have in their favour that be- 
have well when density perturbations come into play: the 
power spectrum for viscous matter is well behaved and 
consistent with large scale structure data [2^]. In par- 
ticular, it does not suffer from the oscillation problem of 
Chaplygin models [2J]. 

On the other hand, since the earliest evidence of tight 
correlations in gamma-ray bursts spectral properties |25| , 
the possibility arose of using GRBs as standard candles. 
Being so GRBs may open a window in redshift as far as 
z r^ 8 |26| . extending then the attainable range provided 
by SNe la observations. For an overview of most recent 
missions, see [27| and references therein. At these epochs, 
z ^ 8, the universe was dominated by dark matter, from 
which follows that this tool is less sensitive to dark en- 
ergy. However for models where dark energy and matter 



are coupled [2^, [2^ or unified Q GRBs might be a useful 
tool. 

GRBs were prevented of being used as standard can- 
dles because the intrinsic faintncss of the nearby events, a 
fact that introduced a bias towards low redshifts of GRB 
and therefore the extrapolation of their correlations to 
low-z events faced serious problems. However this prob- 
lem is solved if one uses SNe la data, in a combined cali- 
bration with GRBs at low redshifts, allowing then GRBs 
be considered as distance indicators (30l - l33j . In particu- 
lar, Liang |34| implemented a method to circumvent this 
objection, obtaining the distance modulus of GRB at low 
rcdshift by interpolating from the Hubble diagram of SNe 
la. Then GRB relations are calibrated without assum- 
ming a particular cosmological model; fair is to say that 
there is some criticism to the method [Sa] . 

In this paper we test the model (presented in [IJ] ) that 
considers a bulk viscous fluid as source of matter, with 
GRBs data given in [36| containing redshifts up to z = 
8.1. In this way we extend the previous SNe la probes 
of the model. In Sec. II we introduced the main aspects 
of the bulk viscosity model; in Sec. HI the calibration of 
GRBs is addressed. In Sec. IV we apply the calibrated 
relation to higher redshifts and do the statistics to find 
the best fit for the bulk viscosity parameter. We carry on 
the analysis with the observational data including SNe la; 
BAG for clusters with redshifts up approximately 2, as 
well as CMB5-year and CMB7-year data from WMAP; 
this last probe tests the model up to the last scattering 
surface, about z ~ 1091.3; according to our analysis the 
model is not reliable to such far epoch, but we extend 
its applicability to epochs of formation of structure z ^ 
8. We discuss conclusions and perspectives in the last 
section. 



II. BULK VISCOSITY DRIVING THE 
ACCELERATED EXPANSION. 

In [13| . the baryon and dark matter components are 
modeled by a pressureless fiuid characterized by a con- 
stant bulk viscosity C. Lying the bulk viscosity constant 
in the range < C < 3, the model possesses many desir- 
able properties, namely: holding the second law of ther- 
modynamics, the derived age of the universe is in perfect 
agreement with the constraint of globular clusters, it has 
a Big-Bang followed by a decelerated expansion with a 
smooth transition to an accelerated expansion epoch in 
late times. Previous cosmological probes done by Avelino 
and Nucamendi include the SNe la Gold 2006 sample 
(182 SNe la) and BAG in a quick fit [13 ■ The many 
good properties of the model justifies to complete with 
GRBs the observational confrontation. In this paper the 
model is tested for the first time using GRBs, and up- 
dated tests using GMB five and seven years WMAP data, 
as well as SNe la Union2 sample and BAG. The results 
point out the strengths and what maybe be the limit of 
applicability of the model, as are the far far epochs of 



the universe, at very high redshifts of the recombination 
epoch. 

Before going on, we note that Nucamendi- Avelino [38| 
addressed the possibility of extending the model of con- 
stant bulk viscosity to a model where the viscosity de- 
pends on iJ as C = Co+Ci^j where Co and C,\ are constants 
and H is the Hubble parameter. Using the SGP Union2 
SNe la data set composed of 557 type la SNe [3^, they 
found that from all possible scenarios predicted by the 
model according to different values of the dimensionless 
bulk viscous coefficients Co and Ci the preferred ones are 
two (Co > 0, Ci < 0) with Co + Ci < 3, and (0 < Co < 3, 
Ci = 0). Nevertheless, one of the disadvantages of the 
model parameterized by C = Co + Ci-^ j is that the best 
estimated total bulk viscosity function C(^) is positive 
for redshifts z <\ and negative for z > 1 latter implies 
a violation to the local second law of thermodynamics 
(see [3^ for details). Then, the simplest model of one 
parameter C = Co turned out to be the best candidate to 
explain the present accelerated expansion. 

For an imperfect fluid, the energy-momentum ten- 
sor with a flrst-order deviation from the thermodynamic 
equilibrium is: 



with 



Tii.v = puy,u^ + {g^,v + UilUv)P*, 



P* =P - CV'u^ 



(1) 



(2) 



where u^ is the four- velocity vector of an observer who 
measures the effective pressure P* ; P and p are the pres- 
sure and density of the fluid respectively. The term C is 
the bulk viscous coefficient that arises in the fluid which 
is out of local thermodynamic equilibrium. The conser- 
vation equation for the viscous fluid is 



u''V^p+{p + P*)V''u^^Q. 



(3) 



In [13 it has been considered a pressureless fluid {P = 
0), then P* — —Cy^u^. We shall consider a spatially flat 
geometry for the Friedmann-Robertson- Walker (FRW) 
cosmology 



ds' 



-dt^ + a^{t){dr^ +r^dn^), 



(4) 



where the function a{t) is the scale factor. The conser- 
vation equation [3] now becomes. 



Pm + (Pm - 3HC)3H = 0, 



(5) 



where H = d/a is the Hubble parameter, pm is the to- 
tal matter density (the dot means time derivative) and 
Vm^ = 3H. 

From the Friedmann equations, the Hubble parameter 
is given through 



Pn^-^HC^Pm-HoC{^) 



1/2 



(6) 



where Hq is the Hubble constant today and it has been 
defined the dimensionless bulk viscous coefficient ^, the 
matter density parameter flmo a-nd the critical density 
today p°,it as 



- ^ 24:TtG ^ ^ PmO 



_ 3ff^ 



Ho 



A^crit 
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so, the Hubble parameter becomes 



ni('.~^-Ui + zf' 



')riO 



(7) 



(8) 



In this model the bulk viscous matter is the only com- 
ponent of the universe implying that the first Friedmann 
equation for a flat universe, H^ = 87rGpm/3, evaluated 
today gives ilmo = 1. With this, the Hubble parameter 
becomes 



H{z) 



Ho r? 



C+ 3-C (1 + ^) 



n3/2 



(9) 



A. Deceleration parameter q 

The function of the deceleration parameter q is defined 



as 



da a 1 

q[a) = -— = — —J. 



(10) 



With the second Friedmann equation, a/a = 
— (47rG/3) {p + 3p), the term a /a can be calculated. For a 
matter-dominated universe with bulk viscosity, it is given 
by 



AttG 



{pm, - 9CH) . 



(11) 



With the first Friedmann equation, H^ = 87rGpm/3, 
we have 



SttG 



H^a). 



(12) 



Then, substituting Eq. [12] and the definition of ( infTTj 
we obtain 



l = l{CHo-Hia))Hi 



a 1 
a 

and substituting [13] into [TO] we have 



^^'■i^'li'-^m 



With 



H{a,C) = 



HofCa^/'+3~c' 



,3/2 



(13) 



(14) 



(15) 



the Eq. [M] is given as 



^(a.C) = 2 



3 - C(l + 2a3/2) 



3-C(l-a3/2) 
then the deceleration parameter today, is given by 



go = g(a= 1,0 = 



1-C 



(16) 



(17) 



Assuming the best estimated values for C, from Table 
[Jwith GRBs, the resulting deceleration parameter today 
is qo = -0.4695 ±0.0324. 

We shall analyze the model presented above. Note 
that this is not the only proposal in this stream, but 
many other similar models exist in the literature, based 
in similar ideas [3, |40| - 149| ; the conclusions we obtained 
then should be pertinent for models alike. 



III. CALIBRATING GRBS 

The main observables that can be measured when 
studying GRBs are its spherical equivalent energy, its 
peak isotropic luminosity, the peak energy of its spec- 
trum, the photon fluence, the energy fluence, the pulse 
duration and the redshift of its host galaxy. Several em- 
pirical correlations among these variables can be estab- 
lished, see for example |50l |. 

The origin, reliability and dispersion of spectral energy 
correlations of GRBs have been debated these late years. 
Some of the pros and cons we comment next on the two 
most used correlations, Ep — Eiso (Amati) relating the 
rest frame energy of the spectra Ep and the isotropic 
energy emitted Eiso, as well as the Ej — Ep (Ghirlanda), 
with Ej = Eisoi^ — COS0 jet) that takes into account the 
non-isotropic release of energy of the GRB. 

Butler et al. [5ll | presented a very critical analysis on 
the observational selection effects, as the origin of the 
GRBs energy correlations. They concluded that rigor- 
ously, neither Ep — Eiso nor Ej — Ep satisfy what should 
a physically-based correlation fulfil: to show a reduced 
scatter in the rest frame relative to the observer frame, 
and that must not persist if the assumed redshifts are 
scattered. Moreover, the discovery of some outliers to 
these relations has raised the suspicion that these cor- 
relations belong only to a sub-population of long GRBs 
or that they are an artifact of the GRBs detection pro- 
cess. In [53 it is pointed out that GRBs correlations are 
strongly influenced by the number-biased against hard 
photons of the GRBs detectors. 

As far as the Ghirlanda correlation Ej — Ep, that in 



some analysis it turns out favored respect to the Amati's, 
as the tighest of the GRBs calibrations, it has the incon- 
venience that to be included in this relation, the GRB 
afterglow must have an observed jet break in its light 
curve, thus, only a fraction of the observed events can 
contribute to establishing this relation. 



Amati [53| contributes to robust the Ep — Eiso relation 
by testing two extremely energetic GRBs (GRB080916C 
and GRB090323) measured by Fermi, showing that both 
events are fully consistent with Ep — Eiso correlation. 
Recently, in [54| it is performed a rigorous analysis of 
the multivariate data, considering also selection effects 
of GRBs, then the authors came to the conclusion that 
there exists a real, intrinsic correlation between Ep — 
Eiso, but not a narrow log- log and it is strongly detector 
dependent. 

Having in mind the previous warnings, we consider the 
Ep — Eiso relation along with data presented in |36| and 
performing the pertinent calibration, we proceed to con- 
strain the bulk viscosity model. 

Most calibrations take for granted a particular cos- 
mological model, to remedy this circularity problem, 
Liang [3J] proposed a cosmology-independent calibration 
method, consisting in calibrating at low redshifts using 
the SNc la data, mending so the few available low red- 
shift GRBs data; the basic idea supporting the method is 
that light travels in the expanding universe in the same 
way no matter what its source be. The typical spectrum 
of the prompt emission of GRBs can be expressed as 
exponentially connected broken power-law, the so called 
Band function [55| . Then we can determine spectral peak 
energy Ep, corresponding to the photon energy at max- 
imum in vFiy spectra. We shall apply the empirical re- 
lation Ep — Eiso that connects Ep 
the isotropic equivalent energy Ei. 



£'p,o6s(l + z) with 
derived by Amati 



Ei. 



47rd^S'6o/o(l 



z)-\ (18) 



where Stoio is the bolometric fiuence of gamma rays in 
the GRB at rcdshift z and cIl is the luminosity distance 
of the GRB (of the hosting galaxy). 

Using at low redshifts the 557 Union2 SNc la data j39| . 
Wei (2010) derived the distance moduli for the 50 low- 
redshift (z < 1.4) GRBs by using a cubic interpolation 
from the 557 Union2 SNe la. Note that with the well- 
known relation 



/i = 5 log 



Mpc 



25, 



(19) 



one can convert the distance modulus fi into luminosity 
distance di, (in units of Mpc) . 

From Eq. [TS] with the corresponding Sboio and lumi- 
nosity distance d^, Wei derived Eiso for these 50 GRBs 
at z < 1.4. Furthermore, with the corresponding Ep for 
these 50 GRBs at z < 1.4, he found the best fit for the 
Amati relation given as 



Wei's calibration was the result of an adjustment called 
bisector of the two ordinary least squares [5g|- Instead 
of using Wei's calibration we preferred to do our own fit 
with a minimum least square method, using the sample 
given in Table 1 by Wei [3g| , consisting of 50 low- rcdshift 
GRBs (z < 1.4). We obtained the following best fit for 
the Amati relation, 

A = 52.7636 ±0.0626, 6 = 1.6283 ± 0.1059, (22) 

The fit is shown in Figure [TJ The errors are calcu- 
lated using error propagation. Our calibration improves 
the one presented by Capozziello and Izzo [57[ , given by 
A = 49.154 ± 0.306, b = 1.444 ± 0.117, performed with 
the SNe la Union sample of 307 objects. Moreover, being 
with errors of better than one part in 10^, our calibra- 
tion is also better than the Amati correlation presented 
by Liang jSJ], done using the 192 SNe la sample. Re- 
mind that at low redshifts the calibration is obtained us- 
ing SNe la samples, and then extrapolated to the GRBs 
with higher redshifts. Therefore the slight differences in 
calibrations are mainly the result of the different SNe la 
samples used at the low rcdshift calibration. 

Extrapolating the calibrated Amati relation to 59 higli- 
redshift GRBs (z > 1.4) given in Table 2 by Wei ^, 
with the corresponding Ep and deriving Eiso from the 
calibrated Amati relation, Eq. [21] (in what follows we 
named MB calibration), the distance moduli fi is ob- 
tained for the extended sample of 59 GRBs at z > 1.4 
using Eos. \rS\ and [TO] and the respective Sboio also re- 
ported in |36| . in order to test the bulk viscosity model. 
To establish a comparison between results obtained using 
Wei's calibration and MB calibration, we include both in 
Tables [1 [ml [V] [VI] [Vll] Otherwise, the GRBs analysis 
was done with the MB calibration. 



IV. DATA ANALYSIS 

The 557 Union2 SNe la data compiled in [33 and the 
59 calibrated GRBs dataset in [3Q| are given in terms of 
the distance modulus ^obs{zi)- The theoretical distance 
modulus is defined by 



( N r, df{z;ai,...,a„) 

^ith(z;ai,...,an) = 5 log — 

Mpc 



25. (23) 



On the other hand, given a parametrization 
H{z;ai, ...,a„) depending on n parameters a,;, we can ob- 
tain the corresponding Hubble free luminosity distance in 
a flat cosmology as 



log- 



Ei, 



erg 



= A + & log 



E„ 



SOOkeV 



(20) 



4''(z;ai,...,a„) =c(l + z) / dz'—— 



Hn 



H{z';ai, ...,a„) 



(24) 



with 



A = 52.7838 ±0.0041 and 5 = 1.7828 ± 0.0072. (21) 



Using the maximum likelihood technique [58| we can 
find the goodness of fit for the corresponding observed 
d'f^^izi). The goodness of fit corresponding to any set 
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Figure 1. 50 GRBs at 2 < 1.4 data in the log (_Bp/300keV) — log (Eiso/erg) plane. The best-fit cahbration is the straight hne 
with slope 6 = 1.6283 ± 0.1059. See Eq. ((22|) . 



of parameters ai, ..., a„ is determined by the probability 
distribution of oi, ..., a„, i.e. 



P(ai,...,a„) = Ne 



r'{ai,...,a,,)/2 



(25) 



where 



X (ai,...,a„) =2^ -2 (26) 



i=l 



(7 



Alois (zi) 



and N IS & normalization factor. The parameters 
oi, ..., On that minimize the x^ expression l26l are the most 
probable parameter values (the 'best fit') and the corre- 
sponding 'X^ifli^ ■■■,an) = Xmin givcs an indication of the 
quality of the fitness for the given parametrization: the 
smaller Xmin i^ the better parametrization. 

For the case of the 557 Union2 SNe la and the 59 cali- 
brated GRBs, x^ is given for the model of bulk viscosity 



xliC,Ho) = 



Uobs(2i) ~ fith{zi,C,Ho) 






M) 



(27) 



Once constructed the x^ function, Eq. 1271 "^6 minimize 
it to find the "best fit" for the free parameters of the 
model. The probability distribution is rewritten as 



P{lHo)^Ne^"'^i^"-^/\ 



(28) 



where N is a normalization factor. 

Besides the SNe la and GRBs, there are other observa- 
tional data very relevant when testing cosmological mod- 
els, so we consider the joint constraints from the latest 
observational data, combined with the SNe la and GRBs, 
namely, the shift parameter R from the WMAP 5-year 
and 7-year data, and the distance parameter A of the 
measurement of the baryon acoustic oscillation (BAO) 
peak in the distribution of SDSS luminous red galaxies. 



In the case of the anisotropics of the cosmic microwave 
background (CMB), the observational parameter that is 
used for the construction of the x? is the shift parameter 
R. Derived in the context of the bulk viscosity model R 
is given by 



RK) = 



E{z',C) ^0 



where the redshift of recombination is z, — 1090.04 
from WMAP5 data [H and z, = 1091.3 from WMAP7 
data [63. 

X^ for i?-CMB is defined as 



XcmbKO 



R{Q — Rob 



(30) 



where i?(C) is the theoretical value predicted by the cos- 
mological model and {Robs, (^R^tJ^ (1.710,0.019) from 
WMAP5 [51 and {Robs, crRjJ)= (1.725,0.018) from 

WMAP7 [eg. 

In the case of baryon acoustic oscillations (BAO), it is 
used the acoustic peak, given by 



A{0 



1 



E^/HZBAO,0 



dz' 



ZBAO Jo E{z'X) 



2/3 



with zbao — 0.35, to construct the function x as 

2 

At/-\ A . 



XBAoiO 



A{C) 



o'A„ 



(31) 



(32) 



where and A{() is the theoretical value predicted by the 
cosmological model. In [6l[, the value of A has been 
determined to be 0.469(n^/0.98)-°-35 ± 0.017. Here the 
scalar spectral index Ug is taken to be 0.963, which has 
been updated from WMAP7 data ^. 



So, the total Xt i^ given by 

Xl ^ Xt{C,Ho) = xl{C,H,)+xhAiBiC)+xlAoiC), 

(33) 
where xl is given in Eq. [571 Xcmb = {R- Robs)^ / o-r and 
Xbao = (^ ^ Aobs)'^/o'A- The best- fit model parameters 
are determined by minimizing the total Xt- 



A. Marginalization over Ho 



Using a prior with the form of a Dirac delta for Hq 
P{Ho) = d{Ho-H*), (37) 

the expression E51 becomes 

P{C) = N ■e-^''-^'"o)/^^ (38) 

where TV is a normalization constant. 



In the statistical process, if some parameters are 
known, this information can be used to 'marginalize' the 
known parameters, i.e. average the probability distri- 
bution [25] around the known value with an appropri- 
ate 'prior' probability distribution. In this work, we 
marginalize over the Hubble constant Hq in order to have 
C as the only free parameter of the model. We use two 
different priors to marginalize Hq: constant and Dirac 
Delta priors. 

1. Constant prior over Hq- 

In this case, we should assume that Hq does not 
have any preferred value a priori, i.e, it has a con- 
stant prior probability distribution function. In- 
stead of minimizing the function x^ given by Eq. 
1271 "wc minimize xip given as 



xlpiC) 



A(C) 



i?(C)+ln(10)/5 



c 



(34) 



where 



A: 



e:Ii (^ 



-" = L»=i — ^ 

r^ — v^n 1 



with 



/i"' 



5 log 



(1 + ^) 



E{z',C) 



25. 



(35) 



(36) 



This new Xcp function docs not depend on Hq any- 
more, the label "cp" stands for constant prior for 
Hq. For more detailed analyzes see the Appendix 
Aof [l|. 

2. Dirac Delta prior over Hq. 

We assume that Hq has a specific value (suggested 
by some other independent observation). In this 
case, the probability distribution has the form of 
a Dirac delta centered at that specific value. In 
particular, we choose Hq = 70.5 km s~^ Mpc~^, as 
reported by 5 year WMAP data, and iJo = 73.8 ± 
2.4km s~^ Mpc~^ as suggested by the observations 
of the Hubble Space Telescope (HST) ^. 



V. DISCUSSION 

In the previous section, we have obtained the con- 
straints of the bulk viscosity model with the latest ob- 
servational data: the combination of 557 Union2 SNe la 
data set in the range 0.02 < z < 1.4 [33], 59 calibrated 
GRBs data set at higher redshifts, the shift parameter 
R from the WMAP5 or WMAP7 data [H H^, and the 
distance parameter A of SDSS luminous red galaxies [6l| . 

As we mentioned above, the best-fit model parameters 
are determined by minimizing the total Xt- For com- 
parison, SNe la and SNe la -I- CMB -t- BAG without 
GRBs have been used to show which is the contribution 
of GRBs to the joint cosmological constraints. In addi- 
tion, some different data sets such as SNe la + BAG and 
SNe la -I- BAG + GRBs have also been jointly consid- 
ered to inspect if the BAG or CMB data are responsible 
of spoiling the statistics. 

In Figure [21 we show the joint confidence regions in the 
(C, Hq) plane for the bulk viscosity model with ( =cte. 
Using 557 SNe la + 59 GRBs (with high redshifts), the 
1-cr confidence region for (C, Hq) of the bulk viscosity 
model is (C, Hq) = (1.9389±0.0647, 69.5616±0.3523 ) 
with Xd.o.f. =0.9572. For comparison, fitting results from 
the data without GRBs are also given in Figure [H (to the 
left). With 557 SNe la, the best-fit values are (C, Hq) 
= (1.9835±0.0668, 69.7130±0.3572) and xl.oj. =0.9830. 
We present in Table [J the best- fit value of ( and Hq with 
1-(T uncertainties, Xmin ^^'^ Xdof^^^ ^^^^ ^^' ^^^ deceler- 
ation parameter qq and the redshift, zt, of the transition 
between the decelerated-accelerated expansion epochs, as 
functions of C; the deceleration parameter is according to 
the ACDM prediction. Moreover, from these values for ( 
and Hq., we see that no significative difference arise from 
probing with SNe la or with SNe la -I- GRBs. Then, the 
validity of the constant bulk viscosity model is extended 
to redshifts of z ~ 8.1. 

In Table [Hi we present the best-fit value of C and Hq 
with l-a uncertainties, and Xmim Xd o f ^ ^^ well as, go 
and Zt but this time using the combination SNe la -t- 
GMB5 -t- BAG, SNe la -I- CMB7 + BAG, SNe la -I- GMB5 
+ BAG + GRBs and SNe la + GMB7 + BAG + GRBs, 
respectively. In all the cases, we found C < that is a 
value unacceptable since it is associated to violation of 
the second law of thermodynamics, however the inclusion 
of GRBs significantly improves the statistics. In Table HIl 




Figure 2. The joint confidence regions in the (-ffo, C) plane for the bulk viscosity model with < (^ < 3. The contours correspond 
to Icr - 4ct confidence regions using Union2 SNe la (left panel) and Union2 SNe la + GRBs (right panel), respectively. GRBs 
have been calibrated using the MB Calibration. The best estimated values and confidence intervals using the Union2 SNe la 
data set are C = 1.9835 ±0.0668 and Hq = 69.7130 ±0.3572 and using Union2 SNe la ± GRBs data set are C, = 1.9389 ±0.0647 
and Hq = 69.5616 ± 0.3523, which are pointed with a dot. 





SNe la 


SNc la ± GRBs | 


Wei's Calibration 


MB Calibration 


c 


1.9835 ± 0.0668 


1.9684 ± 0.0654 


1.9389±0.0647 


Ha 


69.7130 ± 0.3572 


69.6610± 0.3552 


69.5616±0.3523 


Zt 


1.4788 ± 0.1679 


1.4421 ± 0.1573 


1.3726±0.1492 


90 


-0.4918 ± 0.0334 


-0.4842± 0.0327 


-0.4695±0.0324 


Amxn 


544.5880 


570.1470 


586.7870 


2 
Xd.o.f. 


0.9830 


0.9301 


0.9572 



Table I. The best-fit value for the bulk viscosity model parameters (-ffo, C), without priors, with 1— cr uncertainties, Xmm, 
Xd.o./j as well as Zt and go using SNe la, and SNe la ± GRBs, where GRBs have been calibrated using the Wei's Calibration 
or the Montiel-Breton (MB) Calibration. _ffo is in units of km s~^ Mpc~^ and i^, zt and go are dimensionless. The subscript 
"t" stands for "transition" and the subscript "d.o.f". stands for ^^ degrees of freedom" . The confidence intervals are shown in 
Figure [2] 





SNc la + CMB5 + BAO 


SNc la + CMB7 + BAO 


SNc la + CMB5 + BAO + GRBs 


SNo la + CMB7 + BAO + GRBs 


c 


-0.06523±0.03552 


-0.07420±0.03423 


-0.05898±0.03532 


-0.06837±0.03405 


Ho 


62.89780±0.23522 


62.87290±0.23347 


62.88625±0.23501 


62.86002±0.16042 


Zt 


-0.87809±0.04331 


-0.86742±0.03979 


-0.88586±0.04469 


-0.87430±0.04080 


10 


0.53262±0.01776 


0.53710±0.01712 


0.52949±0.01766 


0.53419±0.01703 


Xmin 


1279.08 


1273.84 


1323.06 


1317.91 


2 

Xd.o.f 


2.3005 


2.2911 


2.1513 


2.1429 



Table II. The best-fit value for the bulk viscosity model parameters (-ffo, C), without priors, with 1 

xl.o.f, as well as zt and go using SNe la ± CMB5 + BAO, SNe la ± CMB7 ± BAO, SNe la + CMB5 ± BAO 

SNe la ± CMB7 ± BAO + GRBs, respectively. CMB5 and CMB7 stand for CMB data from 5-years and 7-years from WMAP, 



2 
min 5 

± GRBs and 





SNe la 


SNe la ± GRBs | 


Wei's Calibration 


MB Calibration 


c 

Zt 

go 

2 
Xmin 

Xd.o.f 


1.9838 ± 0.0668 

1.4795 ± 0.1642 

-0.4919± 0.0334 

561.3220 

1.0114 


1.9686 ± 0.0655 

1.4425 ± 0.1576 

-0.4843± 0.0328 

586.8830 

0.9558 


1.9844±0.0666 

1.4810±0.1640 

-0.4922±0.0333 

613.6350 

0.9994 



Table III. The best-fit value of ^ with 1-a uncertainties, and Xmim Xd.o.f. j as well as Zt and go for the model of bulk viscosity 
assuming a constant prior over Hq. These results were obtained using SNe la and SNe la ± GRBs, respectively. 





SNc la + CMB5 + BAO 


SNc la + CMB7 + BAO 


SNo la + CMB5 + BAO + GRBa 


SNc la + CMB7 + BAO + GRBs 


c 


-0.06518±0.03552 


-0.07415±0.03423 


-0.05893±0.03532 


-0.06833±0.03406 


Zt 


-0.87816±0.04332 


-0.86748±0.03979 


-0.88592±0.04470 


-0.87435±0.04082 


90 


0.53259±0.01776 


0.53708±0.01712 


0.52947±0.01766 


0.53417±0.01703 


Xmin 


1296.02 


1290.78 


1340.00 


1334.85 


2 

Xd.o.f 


2.3268 


2.3174 


2.1753 


2.1670 



Table IV. The best-fit value of C, witli 1-cr uncertainties, and Xmim xl.o./.) a-s well as Zt and go for tire model of bulk viscosity 
assuming a constant prior over T^o with SNe la + CMB5 + BAO, SNe la + CMB7 + BAO, SNc la + CMB5 + BAO + GRBs 
and SNe la + CMB7 + BAO + GRBs, respectively. 





SNe la + BAO 


SNc la + BAO + GRBs | 


Wei's Cahbration 


MB Calibration 


c 


1.8693 ± 0.0625 


1.8606 ± 0.0616 


1.8388±0.0654 


Zt 


1.2194 ± 0.1313 


1.2013 ± 0.1279 


1.1566±0.1321 


go 


-0.4347± 0.0313 


-0.4303± 0.0308 


-0.4194±0.0327 


Xmin 


578.138 


603.065 


618.3780 


2 
Xd.o.f 


1.0398 


0.9806 


1.0055 



Table V. The best-fit value of (^ with 1-a uncertainties, and Xmin, Xd.o.f., a-s well as zt and go for the model of bulk viscosity 
assuming a constant prior over Ha using SNe la -|- BAO and SNe la -I- BAO -I- GRBs, respectively. 





SNe la 


SNe la + GRBs | 


Wei's Calibration 


MB Calibration 


c 

Zt 

go 

y^ ■ 

Amxn 

2 
Xd.o.f 


2.0910±0.0429 
1.7662±0.1249 
-0.5455±0.0215 

549.381 

0.9899 


2.0807 ± 0.0423 
1.7365 ± 0.1210 
-0.5404± 0.0212 

575.706 

0.9376 


2.0630±0.0418 
1.6865±0.1162 
-0.5315±0.0212 

593.857 

0.9672 



Table VI. The best-fit value of (^ with l-a uncertainties, and Xmin, Xd.o.f., as well as zt and go were estimated using SNe la 
and SNe la -|- GRBs data set, respectively, for the model of bulk viscosity assuming a Dirac delta prior over Hq located at 
Ho = 70.5 km s~^ Mpc^\ as reported by 5 year WMAP data. 





SNe la 


SNe la -1- GRBs | 


Wei's Calibration 


MB Calibration 


c 

Zt 

go 
v^ ■ 

Amxn 

2 
Xd.o.f 


2.4909±0.0383 
3.5750±0.2764 
-0.7455±0.0192 

671.682 

1.2102 


2.4696 ± 0.0375 
3.4262 ± 0.2534 
-0.7348± 0.0188 

703.588 

1.1459 


2.4440±0.0370 
3.2596±0.2320 
-0.7220±0.0185 

729.296 

1.1878 



Table VII. The best-fit value of C, with l-a uncertainties, and Xmin, Xd.o.f., as well as zt and go were estimated using SNe la 
and SNe la + GRBs data set, respectively, for the model of bulk viscosity assuming a Dirac delta prior over Ho located at 
Ho = 73.8 km s^^ Mpc~^, as suggested by the observations of the Hubble Space Telescope (HST). 





SNc Ia+CMB5-|-BAO 


SNc Ia-|-CMB7-|-BAO 


SNc la -1- CMB5 -|- BAO -|- GRBs 


SNc la + CMB7 -|- BAO -|- GRBs 


c 


2.99393±0.00020 


2.99402±0.00011 


2.99393±0.00013 


2.99402±0.00011 


Zt 


98.09581±2.18115 


99.08961±1. 22986 


98.09581±1.41775 


99.08961±1.22986 


go 


-0.99697±0.00010 


-0.99701±0.00006 


-0.99697±0.00007 


-0.99701±0.00006 


Xmin 


1883.15 


1883.27 


2029.31 


2029.47 


Xd.o.f 


3.3809 


3.3811 


3.2943 


3.2946 



Table VIII. The best-fit value of C, with 1-cr uncertainties, and Xmmi xl o / > as well as zt and go were estimated using SNe 
la + CMB5 + BAO, SNe la -f CMB7 4- BAO, SNe la + CMB5 + BAO' + GRBs and SNe la + CMB7 + BAO + GRBs, 
respectively, for the model of bulk viscosity assuming a Dirac delta prior over Ho located at Ho = 70.5 km s~^ Mpc~^, as 
reported by 5 year WMAP data. 





SNc Ia+CMB5 + BAO 


SNc Ia+CMB7-|-BAO 


SNo la + CMB5 + BAO + GRBs 


SNc la + CMB7 + BAO + GRBs 


c 


2.99439±0.00010 


2.99447±0.00010 


2.99439±0.00005 


2.99447±0.00011 


Zt 


103.45206±1.24359 


104.45889±1.27371 


103.45206±0.62179 


104.45889±1.40107 


go 


-0.99720±0.00005 


-0.99724±0.00005 


-0.99720±0.00003 


-0.99724±0.00005 


Armn 


1641.81 


1641.88 


1775.06 


1775.16 


2 

Xd.o.f 


2.9476 


2.9477 


2.8816 


2.8818 



Table IX. The best-fit value of (" with 1-cr uncertainties, and Xmim xl o / i ^^ ^^ll as zt and go were estimated using SNe la + 
CMB5 + BAO, SNe la + CMB7 + BAO, SNe la + CMB5 + BAO + GRBs and SNe la + CMB7 + BAO + GRBs, respectively, 
for the model of bulk viscosity assuming a Dirac delta prior over Ho located at Ho = 73.8 km s~'^ Mpc^^, as suggested by the 
observations of the HST. 



check first and third column and second and fourth, in- 
cluding GRBs we obtained a better value in Xd o f "The 
negative values of ( has as consequence that go > (see 



Eq. [ni andzt = (2C/(3-C) 
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1 becomes negative. 



Moreover, being C < the pressure turns out to be pos- 
itive since P* = —3(^H. The meaning of this is that we 
can not extrapolate the viscous fluid behaviour to so far 
epochs. 

We can conclude that the model has limited applica- 
bility when extrapolated to very large redshifts like that 
of the recombination epoch, i.e. when the statistical tests 
include the CMB. Notice that all the parameter suffer a 
drastic change in their estimated values, the unreliabil- 
ity of them is measured by the increasing in Xd o f ''^^'^ 
Xmin 3.S 2 and as doubled data, respectively. However, 
the statistics is better when included GRBs than just 
SNe la, i.e. SNe la + CMB + BAO is worst than SNe la 
+ CMB + BAO + GRBs. 

In the following Tables lllli IIVI and |V] are presented 
the best- fit values of C with l-cr uncertainties, and Xmim 
Xd.o.f.' ^■s well as Zf and qo for the model of bulk vis- 
cosity assuming a constant prior over Hq- In Table ITTTl 
are shown the results obtained using SNc la and SNc la 
+ GRBs, which are very similar to those obtained with- 
out prior over Hq. The effect of assuming a constant 
prior on Hq, with and without GRBs, is a best value of 
xlof In Table HV] the results are using SNe la + CMB5 
+ BAO, SNe la + CMB7 + BAO, SNe la + CMB5 + 
BAO + GRBs and SNe la + CMB7 + BAO + GRBs. 
The obtained values are almost the same with or without 
assuming a constant prior on Hq but the best-fit is ob- 
tained without prior, and statistics is better when GRBs 
are included. 

Until now we have seen that when we consider only 
SNe la -t- GRBs we obtain ( > and a reasonable value 
for the Hubble constant Hq. However, when the same 
analysis is performed using SNe la, CMB5 or CMB7, and 
BAO with or without GRBs , we obtain negative values 
of C that disagree with the second law of thermodynamics 
and moreover the estimation for qo and zt is not good and 
also Xdo f g*^*^ worst. 

Since the redshift of CMB corresponds to earlier 
epochs of the universe while BAO refers to more recent 
ones (z ~ 0.35), our guess is that it is the inclusion of 
CMB that ruins our statistics. The aim of Table |V] is 



to prove it. We derived the constraints for the model of 
bulk viscosity using SNe la -I- BAO and SNe la + BAO 
+ GRBs and we found, in both cases, that the value of 
C is in perfect agreement with the value predicted by the 
model to have an accelerated expansion epoch, that is to 
say, 1 < C < 3, and also with this value of C we have an 
accelerated universe in the present with a transition be- 
tween decelerated-accelerated expansion at z w 1.2. This 
result agrees with the analysis on perturbation dynam- 
ics and the calculated matter power spectrum done by 
Zimdahl 2J] on the bulk viscosity model, that turned 
out compatible with the data from 2dFGRS and SDSS 
surveys. 

Finally, we analyze and constraint the viability of the 
model assuming a Dirac delta prior over Hq located at 
Ho =70.5 km s"^ Mpc"^ as reported by 5-year WMAP 
data and at Ho =73.8 km s~^ Mpc~^ as suggested by 
the observations of the HST, the results arc shown from 
Table EI] to Table HXl 

Again, the results obtained with Ho =70.5 km s~^ 
Mpc~^ using SNe la and using the joint SNe la -|- GRBs, 
are in agreement with the model and good values are ob- 
tained for Zt and ^q. However when we use Ho =73.8 km 
s^^ Mpc~^ we obtained values for C, in the expected range 
but with a slight increment in the value oi Xd o f This 
is shown in Table FVll and I VlH On the other hand, when 
we used SNe la, CMB data (WMAP5 or WMAP7) and 
BAO with or without GRBs, assuming a Dirac delta prior 
over Ho located at Hq =70.5 km s'^ Mpc'^ or Ho =73.8 
km s~^ Mpc~^, we found values for ( in the range sug- 
gested by the model but with bad statistic, which is re- 

^3, then 

2/3 

Zt = ( 2C/(3 — C) ) ~ 1 turns out to be enormous, and 
qo ~ —1, for these reasons we conclude that the statistics 
for a Delta Dirac prior on Ho is not reliable. It should be 
noted that when we use SNe la + CMB5 + BAO + GRBs 
or SNe la + CMB7 + BAO + GRBs, i.e. when GRBs 
are included, significantly improve the value of Xd o f ^ 
but anyway, we have a worse fit to data each time that 
CMB is considered. See Tables [Villi and [IXl 

Lastly, if we compare the results obtained using Wei's 
Calibration or MB Calibration for GRBs, as shown in 
Tables [H IIIIl FVl IVII and IVIIl we notice that in most cases 
Xdof i^ closer to 1 when using MB calibration, that is 
to say, the latest calibration provide a better fit. 



fleeted in the value of Xd o f Moreover, since C, 
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VI. CONCLUSIONS 



Modeling the dark sector of the universe with a bulk 
viscous fluid with pressure p = —CO, with a constant bulk 
viscosity (, is tested with several cosmological data: SNe 
fa, CMB and BAO jointly with the calibrated sample of 
GRBs. The model is in good agreement with observa- 
tional data from SN la, BAO and GRBs, up to the range 
of 2 = 8.1 with a chi-square statistics of x^ = 1-01. 

Probing with SNe la and GRBs we obtained as the 
present Hubble parameter Hq = 69.56kms~^Mpc^^, the 
deceleration parameter qq = —0.47, both in good agree- 
ment with the accepted values; as the redshift transi- 
tion parameter we obtain zt = 1.37. When the baryon 
acoustic oscillations (BAO) data are included in the data- 
model confrontation, we obtain a viscosity parameter 
C = 1.84, included in the correct range, and improved 
values for zt = 1.16 and x^ = 1.01. The results with 
BAO confirm the analysis performed by Hipolito-Ricaldi 
(2009) |2J] on the matter power spectrum showing the 
compatibility of the bulk viscosity models with the 2dF- 
GRS and SDSS data. 



However the model docs not pass the CMB test, ex- 
hibiting x^ = 2.15 for CMB5 and x^ = 2.14 for CMB7, 
besides a negative bulk viscosity constant, violating so 
the second law of thermodynamics. These results are 
consequence of neglecting a radiative component in the 
universe, that in so far epochs as the recombination, plays 
an important role. Thus, when testing so far epochs a ra- 
diation component should be added to the universe con- 
tent, a component that is conserved independently of the 
rest of matter considered. 

To conclude, constant bulk viscosity matter reproduce 
observational SN la, GRBs and BAO data in good agree- 
ment (x^ = 1.01), being then a reliable model not only 
for the recent state of the universe, but for epochs where 
structure already existed (z = 8.1). It is an open ques- 
tions if the model could account for the CMB data if a 
radiation component is included in the considered matter 
content of the universe. 
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